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Abstract 

We calculate the two-loop QCD anomalous dimension matrix (ADM) (7^)ndr i n the NDR- 
MS scheme for all the flavour-changing four-quark dimension-six operators that are relevant in 
both the Standard Model and its extensions. Both current-current and penguin diagrams are 
included. Some of our NDR-MS results for AF = 1 operators overlap with the previous calcu- 
lations, but several others have never been published before. In the case of AF = 2 operators, 
our results are compatible with the ones obtained by Ciuchini et al. in the Regularization- 
Independent renormalization scheme, but differ from their NDR-MS results. In order to ex- 
plain the difference, we calculate the ADM of AF = 2 operators again, extracting it from the 
ADM of AF = 1 operators. 



1 Introduction 



Renormalization group short-distance QCD effects play an important role in the phenomenology 
of non-leptonic weak transitions of K-, D- and 5-mesons. An essential ingredient in any 
renormalization group analysis is the anomalous dimension matrix (ADM), which describes the 
mixing of the relevant local four-quark operators under renormalization [2|] . 
The operators considered in the present paper have the form 

^r k A ^^r k B ^i, #?r^#f #f r* (1.1) 

where a, j3 are colour indices and T k A B are generic Dirac matrices given explicitly below. The 
subscripts i in are flavour indices. In the case of FCNC transitions with AF = 2, such as 
neutral meson mixing, one has 

= ^ 2 = (1-2) 

Known examples are the operators (sd)v-A(sd)v-A and (bd)v-A(bd)y^A relevant in the Stan- 
dard Model (SM) to K°-K° and B^-B^ mixing, respectively. 

Four-quark operators that occur in the SM calculations of flavour-changing processes do not 
form a complete set of all the dimension-six four-quark operators. Other operators need to be 
considered in many extensions of the SM, e.g. in the Supersymmetric Standard Model (SSM) 
(see e.g. ref. ||). For instance, the SSM and SM predictions for K°-K° and B®-B® mixing 
can have similar precision only if the two-loop ADM for all the AF = 2 operators is known. 

The main purpose of the present paper is a calculation of the two-loop ADM for all the 
dimension-six flavour-changing four-quark operators in the NDR-MS scheme (MS scheme with 
fully anticommuting 75). Our main findings are the NDR-MS anomalous dimensions of the 
operators with Dirac structures (cf. eq. (|1.1| )): 

r^r^ = (l± 75 )®(l± 75 ) and T k A ®T k B = [^(1 ±75)] ® ± 75 )]. (1.3) 

For these operators, our two-loop results differ from the NDR-MS ones of Ciuchini et al. H, 
but are compatible with their Rl-scheme ADM. For all the other operators, no new calculation 
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is actually necessary — all the two-loop results can be extracted from the existing Standard 
Model ones. 

Our paper is organized as follows. In section |2|, we perform a direct calculation of the NDR- 
MS-scheme ADM of AF = 2 operators. This is a relatively straightforward computation, since 
all the methods are already known from similar SM calculations (see e.g. refs. ||-f|). The 
only novelty here is the introduction of evanescent operators that vanish by the Fierz identities. 

In section |||, we compute the NDR-MS ADM for such AF = 1 operators, to which only the 
current-current diagrams are relevant. Some of the AF = 1 results have never been published 
before. The ones that are not new agree with the old SM calculations. The subject of section ^ 
are AF = 1 operators containing one quark-antiquark pair of the same flavour. We identify 
the operators to which the so-called penguin diagrams are relevant, and give the corresponding 
anomalous dimensions. 

In section |5|, we derive the matrix Af necessary for transforming the Wilson coefficients from 
the NDR-MS to the RI scheme (originally called the MOM scheme) that is more useful for 
non-perturbative calculations of hadronic matrix elements 0. 

Section ^| is devoted to performing a consistency check of our AF = 1 and AF = 2 results. 
The current-current ADM of AF = 1 operators is transformed there to such an operator basis, 
in which the AF = 2 results can be easily read off. This calculation serves also as a preparation 
for the comparison with Ciuchini et al. [Q. Comparison with this article and other existing 
literature is the subject of section [7]. We conclude in section |[ 

In appendix A, we list the evanescent operators relevant to the AF = 2 calculation. In 
appendix B, an analogous list for the AF = 1 case is presented. Appendix C contains additional 
evanescent operators that become important only when one wants to derive the AF = 2 results 
from the AF = 1 ones, as in section |(J Appendix D is devoted to recalling and generalizing 
the notion of "Greek projections" . Finally, appendix E contains a list of separate contributions 
from different diagrams to the one- and two-loop ADMs for AF = 1 operators with Dirac 
structures Q1.3Q - 
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2 Direct calculation of the ADM in the AF = 2 case 



For definiteness, we shall consider here operators responsible for the K°-K° mixing. There are 
8 such operators of dimension 6. They can be split into 5 separate sectors, according to the 
chirality of the quark fields they contain. The operators belonging to the first three sectors 
(VLL, LR and SLL) read 



QX hh = 




Q\ R = 




Q\ K = 


(s a P L d a )(s^P R d^), 


Qf L = 


( s - a p L cr)(^p L ^), 


Qf L = 


{s a a^P L d a ){^a^P L d' 3 ) 1 



(2.1) 

where a^ u = ^[j^Jv] and Pl,r — |(1 T 7s)- The operators belonging to the two remaining 
sectors (VRR and SRR) are obtained from Q\ hh and Qf hL by interchanging Pl and Pr. Since 
QCD preserves chirality, there is no mixing between different sectors. Moreover, the ADMs in 
the VRR and SRR sectors are the same as in the VLL and SLL sectors, respectively. In the 
following, we shall consider only the VLL, LR and SLL sectors. 

In dimensional regularization, the four-quark operators from eq. fl2.1p mix at one loop into the 
evanescent operators listed in appendix A. Specifying these evanescent operators is necessary 



to make precise the definition of the NDR-MS scheme in the effective theory [||, ||, |I0, |TT| . An 
important novelty in the present case (when compared to AF = 1 calculations) is the necessity 
of introducing evanescent operators that vanish in 4 dimensions by the Fierz identities. The 
Fierz identities cannot be analytically continued to D dimensions. Therefore, they have to be 
treated in dimensional regularization in the same manner as the identity 

l^lvlp = g^ulp + 9u P lp, - 9p P lu + ie a pv P 7 a 75, (2.2) 
i.e. appropriate evanescent operators have to be introduced. 
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As an example, consider the operators Qf LL and Q^ L . When these operators are inserted 
into one- and two-loop diagrams, the operators 

gSLL = (s a P L £)(sPP L d a ), (2.3) 

gfL = ^ a ^p Ld P)^a^P L d a ) (2.4) 

are generated. In 4 dimensions these operators can be expressed through Qf LL and Q^ h by 
using the Fierz identities 

{PL)ij{Ph)kl = -{PL)il{PL)kj — -(c I muPl)u(o'' J ' 1J Pl)kj , 
2 o 

(v^PLhi^P^ki = -Q(PMP L )k 3 ~ \((T^P L ) u (a^P L ) kv (2.5) 
which give 

Qf L = 4 -\Qf LL + lQl LL , (2.6) 

Qf L = 4 6Q? LL + igf L . (2.7) 

These relations can be used in the calculation of one-loop ADM. In the case of two-loop calcula- 
tions, in the NDR-MS scheme, where Dirac algebra has to be performed in D ^ 4 dimensions, 
these relations have to be generalized to 

<5f L = -^Qi LL + ^Q2 LL + £i SLL , (2.8) 

Z o 

= ggSLL + IgSLL + f7 SLL (2 _ Q) 

Here, Ef LL and i?| LL are the evanescent operators that vanish in 4 dimensions by Fierz identi- 
ties. They are simply defined by ( |2.8j ) and ( |2.9| ) and are given in appendix A. 

The effective Lagrangian can be written separately for each sector. It takes the form 

(-12 tut2 

Ceff = ! ^(V£V td ) 2 Z$'£C i (jM)[Q i + (counterterms),] , (2.10) 

i 

where Z q is the quark wave-function renormalization constant. 

The coefficients Ci(n) satisfy the Renormalization Group Equation (RGE) 

£ C{p) = 7(//) T d{p) (2.11) 
a/i 



governed by the ADM that has the following perturbative expansion: 

'M = ^7<°> + f$7 (1) + O(o»). (2.12) 

The ADM in the MS or MS scheme is found from one- and two-loop counterterms in the 
effective theory, according to the following relations (equivalent to eqs. (4.26)-(4.37) of ref. |5]): 

7^ = 2a 11 , (2.13) 
7 W = 4a 12 -26c. (2.14) 

The matrices a 11 , d 12 and 6 in the above equations parametrize the MS-scheme counterterms 
in eq. (|2T0|) (for D = A — 2e) 

Oil ^ / 1 OQ 1 



(counterterms)j = 



. k k 



+ (two-loop evanescent counterterms) + 0(a s ). (2-15) 

The matrix c is recovered from one-loop matrix elements of the evanescent operators. Let 
us denote by (-Efc)iioop the one-loop K°-K° amplitude with an insertion of some evanescent 
operator E^. The pole part of such an amplitude is proportional to some linear combination of 
tree-level matrix elements of evanescent operators. The remaining part in the limit D — > 4 can 
be expressed by tree- level matrix elements of the physical operators Qi. The finite coefficients 
of these matrix elements define the matrix c as follows: 



(-E'fe)lloop — — - 



^ dkj{Ej) tree + X! e kj{Fj) tree 
j j 



E C ^)tree + O(e). (2.16) 



Here, Fj stand for such evanescent operators that are not necessary as counterterms for the 
one- loop Green functions with insertions of the physical operators Qi. The matrices c and a 12 
depend on the structure of Fj, but x 1 ' does not. 

The matrices 7^°- ) = 2a 11 , 6 and c in each sector are found from the one- loop ds — > sd 
diagrams presented in fig. [I] with insertions of the physical operators Qi, as well as the evanescent 
operators Ek- We calculate only the "annihilation-type" diagrams, i.e. we drop all the diagrams 
where fermion lines connect the incoming and outgoing particles. Dropping such diagrams 



Figure 1: One-loop diagrams 



consistently at the tree level, at one loop and (later) at two loops does not alter the final results 
for the renormalization constants. 

All the one- and two-loop diagrams considered in the present article are calculated using two 
different methods. In both of them, a covariant gauge-fixing term 

C 9f = -^G;WG:) (2.17) 

is used, and the physical masses are set to zero. In the first method, the external quarks are 
assumed to have momentum ±p. In the second method, the external momenta are set to zero, 
but a common mass parameter is introduced in all the propagator denominators as IR regulator 



12|1 . The two methods give the same results for the MS renormalization constants. The ADMs 
calculated from these renormalization constants with the help of eqs. (|2.13|) and Q2.14Q are 
independent of the gauge-fixing parameter A. 

We begin with presenting the ADM in the SLL sector, because in this very sector our results 
are going to differ (at two loops) from those of ref. [[|. The matrices 7(°) SLL and b SLL are found 
to be the following: 



-6N + 6 + 



6 11 



^(O)SLL UJ ' TU 7 ^ 2 N \ , 218 s 

{ _ 24 -f 2iV + 6-! J ' 



rsLL ( A 



(2.19) 



\ 2N 2 / 

where N stands for the number of colours. 

In order to find the matrix a 12 , we need to calculate two- loop diagrams obtained from the 
ones in fig. [I] by including one-loop corrections on the gluon lines or adding another gluon that 
couples to the open quark lines. Of course, one-loop diagrams with counterterm insertions 



need to be included, too. All the two-loop diagrams and the corresponding colour factors are 
the same as in fig. 2 and table 2 of ref. [f|. However, in the present article, we also consider 



additional Dirac structures ( |1.3|) in the four-quark vertices. 

Inserting the calculated matrix a 12 into eq. ( p. 14 ), we obtain the two-loop ADM. Its entries 
are found to be the following: 

(l)SLL _2roy\r2,lOTy\r,136_12_ 107 , 10 at f _ 2 f _ w_ f 

111 - 6 V_r 3 '3 N 2W 1 3 J 3 J 3N-> ' 

(l)SLL _ _J_AT_ 31 , 9. _ 4 _ J_r , J_ f 
'12 36 iv 9 JV JV^ 18 •> ^9N-li (0 Of]) 

(l)SLL 364 at 704 208 320 _i_ 1_36 r _j_ 1_76 £ 

'21 3 3 N W *~ 3 J 3N-I ' 

7 (1)SLL = SA3 N 2 + 21N _m + U + ^_2^ Nf _ 6f+ _2_ L 

where / stands for the number of active flavours. The above equation is one of the main results 
of the present paper. 

Proceeding analogously in the VLL sector, we reproduce the well-known results for the one- 



and two-loop anomalous dimensions of the operator Q^ LL |J3 



7 (0)VLL _ g _ 



N ' 



(l)VLL _ 19 at _ 22 , 39 _ _57_ , 2 f 2_ £ 

' 6 3 ' N 2N 2 ' 3 J 3N-I ■ 

The matrix b in the VLL sector reads 



(2.21) 



6 VLL =(-5 I). (2.22) 

Finally, our results for the LR sector read 

7 (0)LR = ( * 12 B V (2.23) 



137 _|_ _L5_ _ 22_f 200 at 6 44 f 

~(1)LR / 6 2N 2 3NJ 3 N 3 J \ In nA\ 



>:lr _ [0-5 o 



i 



(2.25) 



As mentioned in the introduction, all the comparisons with existing literature are relegated to 
section [F[ 



S 



3 Current— current contributions to the ADM 
of AF = 1 operators 



In the present section, we evaluate contributions from the current-current diagrams to the 
ADM of AF = 1 operators. For this purpose, we choose the operators in such a manner that 
all the four flavours they contain are different: s, d, u, c. In such a case, the only possible 
diagrams are the current-current ones. 

Twenty linearly independent operators can be built out of four different quark fields. They 
can be split into 8 separate sectors, between which there is no mixing. The operators belonging 
to the first four sectors (VLL, VLR, SLR and SLL) read 

QX hh = {s a i,P L d p ){u^P L c a ) = Q VlVl , 
Q V 2 hh = (s a 7ll P L d a )(u^P L cP) = Q VlVl , 

qVlr = (s« llx P L dP){u^P R c a ) = Q VlVr} 
Q2 LR = {s a ^P L d a ){u?^P RC P) = Q VlVr1 

gSLR = ^ a P L d^P R C a ) = Q LR , 
gSLR = (s a P L d a )(uPP R C?) = Q LR , 

gSLL = (s a P L df)(uPP L C a ) = Q LL , 

QfL = {s a P L d a ){u^P L ^) = Q LL , 

QT L = {s a a^P L d^a^P L c a ) = Q TlTl , 

QT L = {ra^d^o^P^) = Qt l t l j (3.1) 

where on the r.h.s. we have shown the notation of ref. ||. 

The operators belonging to the four remaining sectors (VRR, VRL, SRL and SRR) are 
obtained from the above by interchanging Pl and P R . Obviously, it is sufficient to calculate 
the ADMs only for the VLL, VLR, SLR and SLL sectors. The "mirror" operators in the VRR, 
VRL, SRL and SRR sectors will have exactly the same properties under QCD renormalization. 
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The evanescent operators for the VLL, VLR, SLR and SLL sectors are listed in appendix B. 
Calculation of the renormalization constants and the ADMs proceeds along the same lines as 
in the previous section. The relevant divergences in one- and two-loop diagrams in the cases of 
VLL, VLR and SLR sectors are given in refs. U and JJ. For completeness we give in appendix 
E the corresponding results for the SLL sector. These have not been published so far in the 
NDR-MS scheme. 

Our final results for the AF = 1 ADMs are as follows: 

6 \ 
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(O)VLL 



7 



(l)VLL 



7 



(O)VLR 



7 



(l)VLR 



7 



(O)SLR 



7 



(l)SLR 



7 



(O)SLL 



2V 



6 



2V 



22 _ _57_ 

' 3 22V 2 

« 2V ' nJ 



2 
32V 



/ 



3- 



2V 



-19 N 
6 JV 

22 _ _57_ 

' 3 22V 2 



59 + §/ 



2V 



3 
2 
32V 



203 at2 , 479 , 15 , 10 at f 22 f 
fi JV fi 32V^ 



15 

22V 2 1 3 
_100]V + A + 22f 



-fiV-f +4/ 

137 , _L5_ 22 f 

fi T 9/V 2 3A/J 



_6 

2V 



137 



-6 
-6AT+-J 



18 

A? 



6 ~ 22V 2 
2 JV 

_6_ 

2V 





-f + 24V 
48 



22 f 
3NJ 



3 ^ 2V 



3 +f/ 



+ 4/ 



203 
6 



at2 _|_ 479 j 15 , 

iV + 6 + 22V 2 + 



24 

_48 

2V 



2V 1_ 

2 2V 
1 

2V 4:1 v 



1 

2 

_ J_ 

TV 



fA7 



22 

37V 



./ 



(3.2) 
(3.3) 
(3.4) 
(3.5) 
(3.6) 
(3.7) 

(3.8) 



(l)SLL 2Vj^ I 148 _107_ O AT f ML f 

'11 2 ~r 3 22V 2 ZJV ^ 3JV-'> 

(l)SLL _ 178 )\r , 64 , 16 f 

(l)SLL _ 107 7V 2_71_^_J_ /V /-, J L 

7l3 - 36 ^ V i 8 2V 2 18 iV ^ + 92V' 

(l)SLL _ _ log N s_ _ J_ 

214 - 36 ~ JV 18' 

7 « SLL = -26iV + ^f, 

(l)SLL _ 203 at-2 , 28 107 , 10 Ar f 10 f 

722 - 6" iV + y - 22V 2 + T iV ^ ~ 3JV-/' 



10 



(l)SLL 
/ 23 


18 1 AT 9 J ' 




T (f LL = 
/24 


53 4 _j_ 1 £ 

18 AT 2 ' 9NJ ' 




(l)SLL 
/31 


676 /v 2 _ 1880 — 320 - 

3 3 W 1 


3 J 1 3A^ ^ ' 


7 ^ 2 )SLL = 

/ OZ 


820 at 1 448 88 f 
3 AT 3 J ' 




(l)SLL 

/ oo 


257 A/"2 116 1 21 

18 9 2A" 2 


+ ™ N f + 2f 
' g ,/ ' J ' 


^(l)SLL _ 
734 — 


50 AT - 8 f 
3 iV 3^' 




7ii )SLL = 


488 at , 416 176 f 
3 ' N 3 J ' 




72 )SLL = 


776 320 _|_ 176 

3 A 72 3A^ ^ ' 






3 N 3 J ^ 




7if LL = 


343 at2 i 28 , 21 
18 iV 9 ^ 2A^ 2 


9 1 y J ^ 9N-I ■ 



(3.9) 

Equation (|3.9|) is one of the main results of this work. 



The careful reader has already noticed that the following equalities hold up to o{a s ): 

^VLL _ VLL VLL _ VLL VLR _ SLR 

/ll — /22 J 112 — 121 i 111 ~ 122 i , , 

-,VLR _ „,SLR ,VLR _ ,SLR ,VLR _ ,SLR ^ ' > 

122 — 111 J 7l2 _ 721 > 721 _ 7l2 ' 

At one loop, these equalities are a consequence of the Fierz identities 

(^PlUYPl)m = -h,P L HYPL)k 3 , (3.ii) 

(^PLh(YP R )ki = 2{P R ) u {P L ) kj , (3.12) 

as well as the flavour- and chirality-blind character of QCD interactions. Since the Fierz iden- 
tities are satisfied in four spacetime dimensions only, the relations ( |3.10| ) could be potentially 
broken at two loops in the NDR-MS scheme. Surprisingly, they are not.Q 

On the contrary, analogous relations are broken at two loops in the SLL sector. Because 
of the Fierz relations ( |2.5| ), the one-loop matrix j^ SLL must satisfy the following identity (cf. 
eqs. (9) and (10) of ref. §): 



7 



(0)SLL 



1 In section 4, where the penguin diagrams are considered, no invariance under Fierz rearrangement is 
observed at two loops for the operators with VLL Dirac structure. A detailed discussion of this fact can be 
found in ref. [HI. 
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with 



/ -| |\ 



-I o I o 

6 | 



(3.14) 

V 6 ± 6 / 

No similar relation holds for ^( 1 ) SLL in the NDR-MS scheme. As it has already been said, this 
is not surprising, because the Fierz identities are not true inD^4 dimensions. 

It is unclear to us whether the symmetries (|3.10|) for the VLL, VLR and SLR sectors are 
preserved at two loops in the NDR-MS scheme only by coincidence, or if there is some reason 
beyond this. As we shall see in section || this question is related to the properties of one-loop 
matrix elements of certain evanescent operators. 

4 Penguin contributions to the ADM 
of AF = 1 operators 

In the present section, we shall describe additional contributions to the ADM of AF = 1 
operators that are due to penguin diagrams. Such contributions may arise only when the 
operators contain one quark-antiquark pair of the same flavour. 

For defmiteness, let us consider AS = 1 operators. In the SM analysis of ref. ||, 10 such 
operators were considered 2 



Qi = 


(s a ^P L u^)(u^P L d a ), 


Q 2 = 


(s a 7fl P L u a )(u P YPLd p ), 


Q 3 = 


1 


Qa = 


(s a lfl P L d^J2(Q^PLq a ), 
q 


Q 5 = 


q 


Q 6 = 


(s a lfl p L dP)J2(fYP R q a ), 



2 Our operators here differ from the ones in ref. H by a global normalization factor of 4. Of course, it 
does not affect their ADM. The factor of 4 can be absorbed into the global normalization factor of the effective 



Lagrangian, as the first ratio on the r.h.s. of eq. ( 2.10| ). In this case, the Wilson coefficients of our operators 
are exactly the same as those in ref. 
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Q 7 = -(s a ^P L d a )J2e q (fYP R q P ), 
A q 

1 q 
Q9 = ^(^^Ee^-VPi/), 



Q10 = -(sVLd^J2e q (f^PLq a )- 
A q 



(4.i; 



Their one- and two-loop ADMs, including current-current and penguin diagrams, can be found 
in appendices A and B of ref. ||. They were also obtained in ref. [0. The same results hold for 
the mirror copies of the SM operators, i.e. for the operators obtained from the ones in eq. ( |4.1| ) 
by Pl Pr interchange. 

Beyond SM, new linearly independent operators appear. Their Dirac structures are as in 
eq. (|3.1| ). Our aim is to find a minimal set of linearly independent new operators. In the 
process of identifying these operators, we shall use four-dimensional Dirac algebra, including 
the Fierz relations ( |2.5| ), ( |3.11 ) and (3.12). It turns out that only 3 additional operators (and 
their mirror copies) undergo mixing via penguin diagrams into other four-quark operators in 
eq. (P~T|). These are 

Qu = (sV^ a )[(W) + (A^/ 
Q12 = (s a ^P L d^)[(d^P R d a ) + (s^P R s a j 
Q13 = {s a ltl P L d a )[{d^P R d P ) + {s P YP R s p )\. (4.2) 

The remaining elements of the operator basis can be chosen in such a manner that massless 
penguin diagrams with their insertions vanish. The first three of the remaining operators have 
the structure of Qu, Q13, but with a relative minus sign between the two terms. The next 
two have the structure of Q$ and Qq, but the sum over flavour- conserving currents is replaced by 
a difference between the analogous w-quark and c-quark currents. Their mirror copies have to 
be included, as well. Further operators have the SLL and SRR Dirac structures as in eq. (|1.3|) , 
or they have the form 



[s a P L , R dP)(fP R rf 



(4.3) 
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{s a P LiR d a )(fP R;L qP) 



(4.4) 



where q has flavour different from s and d. It is straightforward to convince oneself that we 
have not missed any linearly independent AS = 1 operator in the above considerations. 

Massless penguin diagrams with insertions of the operators ( |1.3| ), ( O ) and ( |4.4| ) vanish, 
because 

Tr(S odd P L , R ) = and /Y./,S,„:,i /'/../.• = 0, (4.5) 

where S dd is a product of an odd number of Dirac 7-matrices. For dimensional reasons, only 
massless penguin diagrams can cause mixing into other four-quark operators. This means that 
all the AS = 1 operators, except for Q 1 , Q13 and their mirror copies, mix only due to 
current-current diagrams, i.e. their ADMs are identical to the ones we have already calculated 
in sections |2] and 

At the two-loop level, a complication arises because generally the Fierz relations could be 
broken in D 7^ 4 dimensions. Consequently, our use of these relations in the identification of 
linearly independent operators could be put in question. However, as we have already discussed 
in section 2 and will elaborate in section 5, this complication can be avoided by introducing 
appropriate evanescent operators that vanish in four dimensions by Fierz identities. This allows 
us to restrict the basis of new physical operators (undergoing penguin mixing) to the one in 
eq. (f4.2|), even at the two- loop level. 



The introduction of evanescent operators that vanish in four dimensions by Fierz identities 
turns out to have no effect on the two-loop ADM in the case of the operators with VLR and 
SLR structures, because the Fierz identity (|3.12| ) remains valid at two loops in the NDR-MS 
scheme, even if the penguin insertions are considered |J. On the other hand, as pointed out 
in ref. ||, the Fierz identity ( |3.11| ) is broken at two loops in the NDR-MS scheme through 
penguin diagrams, although it remains valid for current-current diagrams. As a result, the 
mixing of the operator 

Q' n = (s a lfl P L d?) [(d^P L d a ) + (s^P L s a )} (4-6) 
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with the operators in eq. pop , through penguin diagrams, differs from the one of Qn at the 
two-loop level. This can be easily verified by using the results of ref. ||. As Q' n = Qn in 
D = 4 dimensions due to the Fierz identity ( |3.11|) , Q' u was not included in the basis ( |4.2|) . By 
working with Qn and the evanescent operator Q' n — Qn, the explicit appearance of Q' u can 
be avoided at any number of loops, so that the basis ( |4.2|) remains unchanged. 

The above discussion implies that the only additional ADMs we need to find in the present 
section are: 



• The 3x3 matrix 7 CC describing the mixing of Qn, Q13 among themselves. 

• The 3x4 matrix 7 p describing the mixing of Qn, Qi 3 into Q 3 , Q 6 via penguin 
diagrams. (Only Q 3 , Q e are generated by massless QCD penguin diagrams with 
four-quark operator insertions.) 

The matrix 7 CC is given by current-current diagrams only. It takes the form 



7c 



7X^=2 

7X^1 



(4.7) 



with 7X^L 2 and *%$L 1 taken from eqs. ( |2T2T| ), Q and Q . 



The matrix 7 P 



7i 0) + 



4vr ,p 

extracted from sections 3.2 and 5.3 of ref. [01. We find 



X 1 ' + ... that originates from penguin diagrams can be easily 



(0) 



7 



Ip 



(I' x ( 1 ' ^' 1 



/ 6iV-| 



352 AT 
27 



16 , 172 

' 3 AT 277V 2 

2 _ 460 

3 27V 

188 



244 1 _20_ 

27 3V 27V 2 



V 



172 at _ 2 1 260 
27 3 ' 27 N 



112 
27 


356 
27V 2 


_32jV 
27 


, 500 
27V 


140 1 
27 ~ t ~ 


148 

27 V 2 


220 N _ 
27 


508 
27V 



22 
3 



6iV + 



3V 



22 
3 



(4.8) 



(4.9) 



The above discussion changes very little in the case of AF = 1 operators, in which F is the 
up-type flavour. Similarly to the AS" = 1 case, all the contributions from penguin diagrams 
can be easily extracted from ref. ||. 
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5 Transformation of the Wilson coefficients 
to the RI scheme 



The ADMs calculated in the present work are given in the NDR-MS scheme that is most 
convenient for perturbative calculations. However, after the Wilson coefficients are evolved 
with the help of RGE (|2.11| ) down to a low energy scale, it might be necessary to transform 
them to another scheme that is more appropriate for non-perturbative calculations of hadronic 
matrix elements |9j . One such scheme is the so-called Regularization-Independent (RI) scheme 
(originally called the MOM scheme) used in ref. ||. Below, we shall give relations between the 
NDR-MS-renormalized and Rl-renormalized Wilson coefficients of all the operators considered 
in sections || and §. 

For completeness, we begin with the definition of the RI scheme. For the massless quark 
propagator, the renormalization condition can be written as 



r 



dp 



sip)-, 1 



= 1, (5.1) 

where /i is the subtraction scale. A simple one-loop calculation is necessary to verify that the 
renormalized inverse propagator in the RI scheme reads 

-.2- 



s ip) R 1 = ~n 



v 



1 _ *L Cf \ (I -In , 
4tt 12 /i 2 



+ 0{a 



2 ) 



(5.2) 



where Cp 



N 2 -l 
2N 



and A is the gauge-fixing parameter (cf. eq. ( |2.17| )). In dimensional regular- 



ization, the corresponding quark wave-function renormalization constant reads 



7R1 



l-^LC F \ (i- 7 + ]n(47r) + i 
47r \e 2 



(5.3) 



provided the subtraction scale \x is identified with the standard MS renormalization scale. 

Conditions similar to eq. (|5.1|) are imposed on renormalized matrix elements of the operators 
( |2.1|) and (|3.1| ) among four external quarks with the same momentum p. The quarks are assumed 
to be massless here. For the AF = 2 operators, such matrix elements have the following form 



(Q\ R )r 



^nV) (Q^tree + if) P»P» ((^7^°) ( ^-fP L dP))^ 

A\fip 2 ) (Q\ R )tr CC + <(P 2 )/P, (isV L d a )is^ u P R ^)) trcc 
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tree j 



/tree 



+ A^( P 2 ) (Q LR ) tree + ^"(^^((^^LfK^WDtree, 

<Q? LL >« = ^ LL (p 2 )(gf L ) trce + ^ LL ( P 2 )(gf L ) tree , 

(gf L ), = AgV) (gf L > trcc + ^ 2 L V) (gf L > tree . (5.4) 

The formfactors Bij(p 2 ) originate from UV- finite parts of Feynman diagrams and are scheme- 
independent. Note that in all the matrix elements multiplied by Bij(p 2 ), only colour-singlet 
quark currents occur. Colour-octet currents are removed from these terms with the help of the 
following Fierz identities (which are independent from the ones in eqs. Q2.5|) , (|3.11|) and fl3.12|) ): 



{pP L ) l3 {pP L ) kl 

(pp L Upp R ) kl 

P^pM W P L )^ VPP B)kl 



(pp L HpP L )kj - -P 2 (i,PLUYP L ) kj , 

2{fP R ) u {pP L ) kj - ~ P 2 (^P R )a(rPL)kr 



No Bij formfactors occur in the SLL sector thanks to the four-dimensional identity 

1 



P»Pu{(?wPl) 1 M UP Pl, 



The RI renormalization condition reads 



^•(-/i 2 ) - u^Bij^-n 2 



with 



)LR 
>21 i 



\ for BY-t h , B\f and B\ 
otherwise. 

The renormalization condition (|5.9|) can be equivalently written as 



A 



effective / 2 



-/x 2 ) = 8ij, 



with ^ffective obtained f r0 m eqs. 

(pP L ) ® (pP LjR ) 
p^(a^P LjR ) ® (a up P R , L ) 



) by making the following ad hoc replacements 

\p 2 (^PL)®(rp L , R ), 

0. 



(5.5) 
(5.6) 
(5.7) 

(5.8) 

(5.9) 

(5.10) 

(5.H) 



(5.12) 
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In the case of AF = 1 operators, the general structure of one-loop matrix elements is similar 
to that in eq. (|5.4j) , but the number of formfactors is larger, because operators with colour-octet 
currents are now linearly independent. The matrix elements can be written as 

{Q t ) R = Af ct ™(p 2 )(Q 3 ) trce + Ni, (5.13) 

where iV, vanish under the replacements (|5.12|) . The RI renormalization condition then has the 
same form as in eq. ( |5.11[) . 

In each of the sectors, the Rl-renormalized Wilson coefficients can be obtained from the 
NDR-MS-renormalized ones with the help of the following relation 

={l- ^AfSg^fo)) C^) + O (a 2 ) , (5.14) 

where 



AtT r i 

[tovBuMv = \ A W) ~ <V)] ■ ( 5 - 15 ) 

The above relations can be easily derived from the fact that the renormalized matrix element 
of the whole effective Hamiltonian is scheme-independent, i.e. 

C m (»)(Q(n,p 2 )r = C m (ri(Q(v, P 2 )) m . (5.16) 

Again, the RI subtraction scale and the standard MS renormalization scale have been tacitly 
identified. The external states must be the same in eq. ( p,16| ). Consequently, the Rl-scheme 
renormalization constant (|5.3| ) must be used for external quark lines in Afj S (p 2 ) that enters 



into eq. ( |5.15| ). 

The dependence on p 2 and the explicit dependence on \i cancels out in Ar^^ RI ( |5.15|) . 



However, one should not forget that this matrix depends on the gauge-fixing parameter A that 
is, in turn, /i-dependent. 

Once the RI renormalization conditions have been specified, finding the explicit form of 
^ms^ri i s om y a matter of a straightforward one-loop computation. Our results for the 
AF = 2 operators follows: 

A vll 7_ X _ i21n2 + 121n2 + A(^-^--41n2 + ^) 

^'MS^RI ' N l-tUlZ^- N -\- A ^ 2N & \ AT y ' 
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^'ms->ri 




Ar^- L 

^' MS^RI 

Ar^- L 

^' MS->RI 



Ar- 



SLL 

MS^RI 



11 



12 



21 



Ar^- L 

^' MS->RI 



22 

In the AF = 1 case, we find 



+ ^¥ + X {m + ^¥) 4 + 41n2 + A(l+4hi2) 

■1 + 1b2-A(|-1h2) -AN + f + ^ - A (f - £ - ^) 

-4iV + 7 + A_41n2 + ^ + A(i + ^-f + ^), 

_5_ 13 _ 2 In 2 , 5 In 2 i \ f _5 J ln_2^ , ln_2\ 

12 12AT 3 6N A \24 6N 3 ~r 6N J ' 

4 _ ^ _ 32 In 2 + - A (2 + f + 16 In 2 - ^) , 

26 In 2 _i_ \ f A, 7 §_ 

3W A ^ 2 6 6W 



1 5 28 In 2 _|_ 26 In 

3 3N 3 + 3JV 



81n2 | 10 In 2 \ 
3 + 3if j' 



^' MS-+RI 



^' MS-»RI 



^' MS-+RI 



Ar^- L 

^' MS->R 



Ar- 



Ar- 



SLL 

ms^r: 

SLL 
MS->R 



A SLL 

^' ms^r: 

ms^r 
Ar^- L 

^' MS^R 



Ar- 



SLL 

MS^R 



MS^R 
A SLL 

^'ms^r: 
Ar— 

MS->R 



Ar- 



Ar- 



SLL 

ms^r: 

SLL 

ms->r 



a SLL 
^' MS-»R 



11 
12 
13 
14 
21 
22 
23 
24 
31 
32 
33 
34 
41 



" N 



7 , 12 In 2 _ \ f_3 41n2\ 

~T N \2N N ) 

- 121n2 + A(| -41n2) 

_4 A T + | F+ 21n2_ X f M 



7 



7- 121n2 + A(| -41n2) 

10 1 ■' ' " 41n2 \ 



12 In 2 _ \ ( _3_ 

N \2N 



N 



I 2 In 2 \ 

iv iv " V2W AT J 

2-21n2 + A(l-21n2) -4iV + | 

2 In 2 , \ fj_ , 21n2 \ 
AT "I" A V2Af N ) i 



21n2 + A(l -21n2) 

2 In 2 I \ / 1 I 21n2 \ 
fv ~ t ~ A V2A r ~ t ~ A 7 / 

+ 21n2) 



A 



_3N _i_ _5_ _i_ 21n_ 

2 AT N ~ * \2N ~ 

= -I -21n2- A (| + 21n2) , 

N_ _ 13 , 5 In 2 , \ ( N 1_ 

2 12AT 6N ~r A V 8 6N 

J 5 In 2 , \ ( \ ln_2\ 

12 6 + \24 6 ) ' 

= -l-21n2 + A(l-21n2), 

^ at I 5 | 21n2 \ / 3N _L 21n2 A 

A 1 " Af A \ 2 2Af A" / ' 



T 6 AT J 



13 
12 



_5_ i 2 In 2 

N N 

5 In 2 , \ ( 1 In 2 \ 



13 , 5 In 2 _ \ ( _L ln2\ 

" 12N 6AT A V 6Af 6AT J ' 

= 4iV-i + ^ + A(6iV-| + ^), 



AT N y 

40 ln2 + A(2-81n2), 

5JV 5_ , 26 In 2 \ f Af \ 5 

2 37V 3N V 6AT 



2 3AT 

25 26 In 2 , \ / 11 
~6 3~ + A [~6 ~ 

12-40 ln2 + A(8-81n2) 



10 In 2 \ 

3N J ' 



10 In 2 \ 
3 J ' 



19 



„SLL 
MS- 

„SLL 



Ar- 

MS-+RI 



Ar- 

' MS^RI 

Ar^- L 

MS-*RI 



12 , 40 In 2 _ \ (S_ 
' N "T" AT A 

+ A(1-^) 

5_ , 26 In 2 i \ f N 5_ , 10 In 2 

I , 3iV 3iV \ 2 6N ' 3iV 



In 2 

42 x x \ N ) 1 
5 _ 26 In 2 

43 3 3 1 ' v V3 3 7' 



In section 0, the above results will be used in performing the comparison with ref . [f| . 

6 Recovering the ADM of AF = 2 operators 
from AF = 1 results 

Let us now use our AF = 1 anomalous dimensions from section § to find again the ADM of 
AF = 2 operators. This will serve as a cross-check of our findings and as a preparation for the 
comparison with ref. M in section [T[ 



Starting from eq. (|3.1| ), we shall pass to another operator basis where the operators are either 
symmetric or antisymmetric under d c interchange. Next, the flavours of both quarks and 
both antiquarks will be set equal. For definiteness, we shall do it first in the SLL sector. The 
superscript "SLL" will be understood for all the relevant quantities below, and we shall not 
write it explicitly. 

In four spacetime dimensions, passing to the new operator basis would be equivalent to 
performing a simple linear transformation of the operators. In the framework of dimensional 
regularization, introducing additional evanescent operators becomes necessary. In the SLL 
sector, only two evanescent operators were needed in the AF = 1 calculation (see appendix B). 
Now, we need to introduce six additional evanescent operators in this sector. They are defined 
in appendix C. 

We begin with a redefinition of the physical operators Qi (i = 1,...,4) that amounts to 
adding to them appropriate linear combinations of the evanescent operators Ec 



Qi^Qi + J2 w ikE k = [Q*Lw (6-1) 
fe=i 



20 



where 

/00-±0±000^ 
00 00000 
00 6 ± 
\0 00000/ 

The "new" operators read 



W 



(6.2) 



[Qi] 
[Qs] 



-\[Qi] F + ^[Q^ F , 
Q[Qi] F + l[Q3} F , 



IQ2) 

[Qa] 



Q2, 



where 



[Q 



(s a P L c a ) (u p P L 



[Q 



3Jf 



(6.3) 
(6.4) 

(6.5) 



In 4 spacetime dimensions, the transformation fl6.1| ) would be equivalent to performing the Fierz 
rearrangement of Qi and Q3, as E}. would not contribute. Since the Fierz identities cannot 
be analytically continued to D dimensions, the Fierz rearrangement must be understood in 
terms of the transformation (|B~T| ), so long as the MS scheme is used. The MS-renormalized 
one-loop matrix elements of Q\ and Q3 are affected by this transformation. This means that 
the renormalization scheme is changed. We pass from one version of the NDR-MS scheme to 
another, even though the evanescent operators remain unchanged. 

After the redefinition ( |6.1| ), we perform a simple linear transformation of the operators 



[Q 



*Jnew 



H Ri i [Qj] 
j = l 



(6.6) 



with 



R 



1 1 

'4 2 



V -3 



_L 

16 
1 
4 

_ J_ 

16 
_ 1 
4 



1 
2 





(6.7) 



As one can easily check, our final operator basis is {Qf, Q£, Qi , Q2 }■> where 

1 

2 



Qt 



{s a P L d a )(u p P L c p ) ± (s a P L c a )(u p P L d p ) 
[s a ^ u P L d a ){u^^P L c^) ± {s a a^P L c a ){u^^P L d p ) 



{6.1 
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The ADM transforms as follows: 



7 



(o) 



7 



(i) 



R { 7 « + [Ar, 7 (0) ] + 2/3 Af } R~ l 



(6.9) 
(6.10) 



where (3q = — |/. The matrix Af reflects in the usual manner [|T4| change of the renor- 



malization scheme that follows from eq. (|6.1| ). The explicit form of Af is 



Af 



-WL 



(6.11) 



provided We = 0. The matrices c and e are found from one-loop matrix elements of evanescent 
operators, as in eq. ( |2.16| ). The product We indeed vanishes in our case, and 



V 



3jV — — 

4 iv N 

77V _ 28 
(iv N 



17 
4 
★ 

21 
★ 



J- TV — 

16 4JV 

4 iv ^ N 
★ 



_3_ 

16 

_ 13 
4 



(6.12) 



Here, stars denote non-vanishing elements of c that are irrelevant for us, since they do not affect 
the matrix 



' 2 iV ^ AT 



Af 



-Wc 



\ 



SN + f 




__ l 

2 



-36 




4 JV 4iV 



I_/V — — 







(6.13) 



/ 



After the transformation (|6.9| , |6.10|) , the ADM in the basis {Qf, Qt, Q\ , Q2 } is found to have 
the form 



74x4 

where 7—7 
7 (o)± 



72x2 2X 2 \ 
02x2 72X2 / ' 

(0)± + f^A(l)± + 

An 



-6N ± 6 + f 
T24-f 



(6.14) 



±1 _ JL 

2 Af 

2iV±6-| 



(6.15) 
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and 



(i)± 
(i)± 

7l2 

(i)± 

721 

(1)± 

722 



203 /V 2 + I07/V-l-iM = ri 2 ._ i07 , 10 at j? _ 2 i- 10/ 

6 3 " 1 " 3 + JV 27V 2 "l" 3 ^ T 3 ^ 3JV-' ' 

^36 9 N N 2 ' 18 J ~ 9N ■> ' 

— t— 364 at 704 -j- 208 320 i 136 r , 176 f 
+ 3 3+W TV 2 - 1 - 3 ■> 3N J ' 



(6.16) 



One can easily verify that the matrix 7 + is equal to the one we have already found in 
eqs. (|2.18| ) and (12.20) ). It must be so, because the operators Qf from eq. fl6.8|) reduce to 
Qf hL from eq. ( |2.1| ) when the flavour replacements c d and u — > s are made. Moreover, the 
evanescent operators listed in appendices B and C can be linearly combined to the ones that are 
either symmetric or antisymmetric under d <-> c interchange. When the flavour replacements 
c — > d and w — > s are made, the antisymmetric operators vanish, while the symmetric ones 
become equal to those in appendix A. Thus, we have shown how to extract the AF = 2 results 
from the AF = 1 ones. 

Let us now briefly describe the analogous transformations in the VLL and LR=VLR©SLR 
sectors. All the necessary evanescent operators are given in appendices B and C. The relevant 
matrices W and R are the following: 



W 



VLL 



R 



VLL 



I 1 
^ 

2 ( -1 1 )' 



It 



LR 

4x12 



vv 2x6 
02x6 



U2x6 
TI/SLR 
vv 2x6 



-2W YLL 

02x6 



U2x6 
-|W VLL 



R 



LR 



/ 1 -2 \ 



-| 1 




1 

2 



1 2 





1 / 



(6.17) 



Consequently, the final operator bases are {Q\ LL+ , Q\ LL ~} and {Qi R+ , Qi R ~ , Q% R ~}, 



where 



Ql 



VLL± 



[s^PLd^YPLC 13 ) ± (s^PLC^YPLd?) 
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Q L ! R± = ^[(^P L d a )(uP^P R cP)±(s a 1 ,P R C a )(uPYP L df 3 ) 

Q L 2 R± = \[{rP L d a ){uPp R cP)±{s a P R c a ){uPp L d?)]. (6.18) 

An important simplification in the present case is that the one-loop matrix elements of the 
evanescent operators E^ LL , E^ LR and i?| LR from appendix C vanish in the limit D — > 4, 
after subtraction of the MS-counterterms proportional to evanescent operators only. This 
means that the third rows of c VLL , c VLR and c SLR vanish (cf. eq. ( |2.16| )). Consequently, 



Af VLL = _jyVLL £ VLL = q and Af LR = — j^/LRgLR = q Thig ^ why ^ twQ . loop 

AF = 1 matrices of the VLL, VLR and SLR sectors exhibited Fierz symmetry in eq. ( |3.10| ). 
The transformations of the two-loop ADMs in the VLL and LR sectors thus look as if we worked 
in 4 dimensions, i.e. they reduce to simple multiplications by the corresponding .R-matrices 
and their inversions. The final results are 

"VLL (l VLL+ \ LR / 72^ + 2x2 

72X2= n A.VLL- ], 74X4= „ . r p_ ) l 6 - 19 ) 



/ V U 2x2 7 2 x2 



where 



i mLL± = ±6 -A 

= ^iv-^i^-^i^-A;, 



6 3 iV 2iV 2 3' 3iV' 
±12 



~(0)LR± __ I N . 

\0 -QN + §)' 

/ 137 , _L5_ _ _22_ r 4-200 A/" xixMf \ 

-(1)LH± _ | 6 ~T 2V 2 3JV-' 31 3 iV ± V ± 3 ■> 1 ffi 9fV> 

One can see that 7 V ' LL+ and 7 iR+ are identical to our AF = 2 results in eqs. ( 2.21|) , ( |2.23|) and 
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7 Comparison with previous ADM calculations 

In the present section, we compare our findings from sections 0, |3] and |] with the previously 
published results for anomalous- dimension matrices. 

7.1 One- loop results 

As far as the one-loop QCD ADMs of four-quark operators are concerned, the historical order 
of their evaluation was as follows: 

• Current-current contributions to the one-loop ADM of AF = 1 operators belonging to 



the VLL and VLR sectors were originally calculated in refs. [15|, [16|]. These results were 
also immediately applicable to the SLR sector, because the Fierz rearrangement has a 
trivial effect at one loop. For the same reason, one-loop anomalous dimensions of the 
AF = 2 operators belonging to the VLL and LR sectors could have been immediately 
read off from these articles. Thus, after 1974, the only unpublished one-loop current- 
current anomalous dimensions were those of the SLL sector, both in the AF = 1 and 
AF = 2 cases. 

One-loop penguin contributions to the ADM of the Standard Model operators were origi- 
nally evaluated in refs. |17]]-[T9|]. As we have shown in section |], penguin contributions to 
the ADM of other (beyond-SM) flavour- changing dimension-six operators can be easily 
extracted from the SM calculations, both at one and at two loops. 



• To our knowledge, the first published results for ^°'> SLL occur in refs. |(| and @, for the 
AF = 2 and AF = 1 cases, respectively. 

The one-loop ADMs given in the present article agree with all the papers quoted above. How- 
ever, in order to perform comparisons, one often needs to make simple linear transformations, 
because different operator bases are used by different authors. For instance, the results for 



'po) in ref. [ 20 1 are given in the basis {Q\ ,Q\ }■ In order to compare them with our 



eq. ( |2.18| ), one should use the relation (p^6|). Similarly, eqs. ( |6.7| ) and ( |6.9| ) need to be used for 
comparing our ^(°) SLL in eq. ( |3.8|) with the corresponding results in ref. ||. 
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7.2 Two-loop results 

The history of previous two-loop computations is as follows: 



• The current-current anomalous dimensions of the AF = 1 operators belonging to the VLL 
sector were originally calculated in ref. (in the D RED-MS scheme), and confirmed in 
ref. [||] (where the NDR-MS and HV-MS results were also given). 

• The remaining elements of the two-loop QCD ADM for AF = 1 operators relevant in the 
SM were calculated in refs. || [7]. New results in these papers were the current-current 
contributions in the VLR sector, as well as all the penguin contributions. The SLR sector 
results in the AF = 1 well as the AF = 2 results for the VLL and LR sectors 
could be easily derived from them with the help of Fierz identities, because the NDR- 
MS-renormalized one-loop matrix elements remain invariant under Fierz transformations, 
except for the current-current ones in the SLL sector, and the penguin ones in the VLL 
sector. Therefore, in the early 1990's, the only unknown two-loop anomalous dimensions 
were those of the SLL sector. 

• The first calculation of the two-loop ADM in the SLL sector was performed by Ciuchini 
et al. Q, in both the AF = 1 and AF = 2 cases. The ADM was calculated there in 
the so-called "FRI" renormalization scheme. The transformation rules were given to the 
LRI scheme (Landau-gauge Rl-scheme) and to the NDR-MS scheme. Current-current 
anomalous dimensions for the remaining sectors were recalculated as well. 

• Penguin contributions to the ADM of non-SM operators are considered for the first time 
in the present article. 

All the two-loop results presented here agree with the previous calculations mentioned above, 
except for the NDR-MS ones for the SLL sector found in ref. [|J. Below, we explain the reason 
for this disagreement. 
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7.3 Comparison with ref. 



In ref. Q, the two- loop ADM for AF = 1 operators of the SLL sector was given in the basis 
defined in eq. (13) of that paper, which is equivalent to our eq. fl6.8|). It was presented in the 
so-called "FRI" renormalization scheme, and the transformation rules to the NDR-MS scheme 
were appended. Applying these transformation rules to their "FRI" -scheme ADM, one obtains 
results that differ from our eq. QS.16| ). In particular, a mixing between and Qf occurs, which 
is absent in our result ( 6.16Q . We could obtain their result if we ignored the transformation 
fl6.1| ) and, consequently, used Af = in our eq. fl6.10| ). However, the final results would then 
correspond to the basis 



Q't 



Q't 



2 o 



[s a *^P L d a )(uP*» u P L cP)±6{s a P L dP){u f3 P L c a ) ± -{s a a, v P L d^a^P L c a ) 



(7.1) 



rather than the one in eq. ( |6.8| ). In 4 spacetime dimensions, the operators ( |6.8|) and (71 
are identical, thanks to the Fierz identities (|2.5| ). However, in D dimensions they are not. 



Consequently, their NDR-MS-renormalized matrix elements differ at one loop, and it is not 
surprising that the two-loop ADM depends on which of the two bases is used. 

We informed the authors of ref. |4[] about our findings prior to publication of the present 
article. They responded that although their NDR-MS results had been claimed to correspond 
to the basis Q5.8|) , the NDR-MS renormalization conditions had been actually imposed in the 
basis (|7.1|) . However, they had forgotten to mention this in their article. Unfortunately, such a 
mistake in the presentation has the same effect on the final result as a mistake in the calculation 
that amounts to missing Af ^ in eq. (|6.10|) . 

As far as the two-loop ADM for AF = 2 operators of the SLL sector is concerned, the 
situation is as follows. If we made the flavour replacements c — > d and u — > s in the basis (|7. 1|) , 
but did not change anything in the ADM, we could interpret this ADM as the one for AF = 2 
operators, as the authors of ref. [|J did. However, it would correspond to quite non-standard 
conventions for the treatment of the evanescent operators obtained from Q'± and Q'2 after 
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the flavour replacements. One would need to assume that the finite one-loop matrix elements 
of these evanescent operators are not renormalized away, contrary to the usual procedure for 
any evanescent operator |5|, [|, [H| |ll[] . Such non-standard conventions make the RGE evolution 
more complicated, because one has to deal with a 4 x 4 instead of a 2 x 2 ADM in the NDR-MS 
RGE for the SLL sector, in the AF = 2 case. The calculation of the one-loop matrix elements 
becomes more involved, as well. 

In the AF = 2 case, no calculation is necessary to convince oneself that the results of 
ref. [|J cannot correspond to the NDR-MS renormalization conditions imposed in the basis 
(|6.8|) (their eq. (13)). Once the c — > d and u — > s replacements have been made, the operators 
in eq. ( |6.8|) vanish identically in D dimensions. Therefore, they cannot mix into the Qf 
operators, independently of what the treatment of evanescent operators is. On the other hand, 
mixing of into Qf was claimed to be found in the NDR-MS scheme in ref. [Q. Therefore, 
an inconsistency is clearly seen. 

In the remainder of this section, we shall verify that our NDR-MS results are compatible 
with the LRI ones of ref. [|J. By differentiating eq. (|5.14| ) with respect to fi, one obtains 

f3 a 2 s (ii) 



in 



Af^— 

2 MS^RI 



+ 1 



(7.2) 



where we have used the RGE (|2.11| ) 



H—ol s (ijl) 
djj, 



2tt 



+0(al 



and 



2tt 



(7.3) 



We have also used the fact that the dependence of Afjj^^ RI on \i originates solely from its 
dependence on the gauge-fixing parameter A(/i). 

Next, we use eq. ( |5.14j ) again to express C MS (fi) by C RI (/i) in eq. Q7T2|). Then, the first two 
terms of the perturbative expansion (|2.12|) of 7rj can be easily read off 



off 
off 



1- 



MS' 



7^ + 

'MS 



(7.4) 
(7.5) 
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Armed with our explicit expressions for Ar^^ RI given in section |5] and with the values of 

= and (7.6) 

we can easily calculate the Rl-scheme ADM from our MS results, for arbitrary A. Setting then 
A — ► 0, we recover all the LRI-scheme anomalous dimensions given in ref. Q. 

As far as the "FRI" -scheme ADMs of ref. [[| are concerned, we can confirm them as well. 
However, it should be emphasized that the "FRI" scheme is not equivalent to the RI scheme 
considered in section [| for any choice of A. The "FRI" scheme cannot be defined beyond 
perturbation theory, because different external momenta are chosen in different diagrams when 
the renormalization conditions are specified. Therefore, in our opinion, the main advantage of 
the RI scheme is lost. 



8 Conclusions 

In the present paper, we have calculated the two-loop QCD anomalous dimensions matrix 
(ADM) (j( 1 ')ndr in the NDR-MS scheme for all the four-fermion dimension-six flavour- 
changing operators that are relevant to both the Standard Model and its extensions. 

The AF = 2 two-loop results can be found in eqs. ( ggOfl , ( ggg ) and ( ggg ). While the 
matrices in eqs. ( 2.21Q and ( 2.24 ) could be extracted from the already published results, the 



two-loop NDR-MS ADM ( 2.20 ) for the SLL operators defined in eq. Q2.1|) is correctly calculated 
for the first time here. 

The AF = 1 two-loop results for operators containing four different quark flavours can be 
found in eqs. (p]), (|3j), Q and (p|). While the matrices in eqs. Q, fl3j) and Q could 
be extracted from the already published results, the two-loop NDR-MS ADM ( |3.9| ) for the SLL 
operators defined in eq. (|3.1|) is correctly calculated for the first time here. 

Penguin contributions to the ADM of non-SM operators have been considered for the first 
time here. These contributions can be easily extracted from the existing SM calculations. 
We have identified the relevant non-SM operators in the AS" = 1 case, and presented the 
corresponding ADM explicitly in eqs. ( |4.8|) and 
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We have demonstrated that the main findings of our paper, given in eqs. ( |2.20| ) and ( |3.9| ). 
are compatible with each other, i.e. we have shown how to properly transform the ADMs from 
the AF = 1 to the AF = 2 case. In this context, we have pointed out that in the process of 
this transformation it is necessary to introduce additional evanescent operators that vanish in 
four spacetime dimensions because of the Fierz identities. 

We have also given the rules that allow transforming our NDR-MS ADMs to the correspond- 
ing results in the RI scheme, for arbitrary gauge-fixing parameter A. They can be found in the 
end of section 5. 

The AF = 1 two-loop ADMs for all the operators defined in eq. fl3.1|) were previously 
presented in ref. |4j], in the Qf basis. In the case of VLL, VLR and SLR operators, there is full 
agreement between their and our results. The case of SLL operators is more subtle. We can 
confirm their LRI-scheme results (RI scheme with A = 0). However, their NDR-MS ADM is 
compatible with ours only after correcting their eq. (13), i.e. after changing the definitions of 



their SLL operators to the ones given in eq. (JO 



After such a correction in eq. (13) of ref. also their AF = 2 NDR-MS results are 
compatible with ours, provided they are understood in terms of quite non-standard conventions 
for the treatment of evanescent operators. In their conventions, the two-loop AF = 2 NDR-MS 
ADM is a 4x4 rather than 2x2 matrix, which makes the RGE evolution and calculating low- 
energy matrix elements unnecessarily complicated. Consequently, the results presented here 
should be more useful for phenomenological applications. 
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Appendix A 

Here, we specify the evanescent operators that are necessary as counterterms for one-loop 
diagrams with insertions of the AF = 2 operators ( |2.1|) . 



77. VLL 

n, x — 


('« a o/ Pr/VeV P T rl a \ n VLL 


771 VLL 


is 7/i7w7p-ri« Jl s 7 7 7 "i« J + lo + 4ej(^ 1 , 


77. VLL 
^3 


fe a -y <-y Pr tf^ V S^VVV'Pr (f*") + f - 1 6 4- 4fl OY LL 
V" 5 //x /V /p- 4 L u A" 5 / / / ^ u / ' v ~ ^ C /Vl ) 


E\ R = 


(s a P L £)(s p P R d a ) + -Qtj R , 


£ 2 LR = 


(s a ^aPLd p )(s p ^P R d a ) + 2Q^ R , 


p LR _ 


1-5 7 /i 7;/7px7 / ct J(s 7 7 7 -Trcs j + ^—4 — 4ej(^ 1 , 


i? 4 LR = 


(s q 7m7,7p^)( A M 7VV) + (8 + 8e)g^ R , 


£ 5 LR = 


(sV^P^X^P^) - 6eQ^ R , 


^ = 






(s a P L d^P L d a ) + l-Ql LL - ]-Ql LL , 


£<SLL _ 


{s a a^P L d^a^P L d a ) - 6Qf L - ^ LL , 


^rSLL _ 


(s%7*7p7^ a )(AV7V^) + ("64 + 96e)Qf l 


^rSLL _ 


(s a W»JpJ* p Ld P )(^YYj CT PLd a ) - 64Qf L + (-16 4 



-16 + 8e)Q^ 



SLL 



iSLL 



The evanescent operators for the VRR and SRR sectors, i.e. and E^ RR are obtained 

by replacing L by R in the definitions of E^ LL and Pf LL . 

The operators E^ LL , P]" R , E^ R , Ef LL and Pf LL vanish in four spacetime dimensions because 
of the Fierz identities (ETTlD , (|]T|) and The operators P^ LL , £^ LL , E^ R , E\ K , Pf LL 

and Pf LL vanish by the four-dimensional identity (|2.2|) . Finally, E^ R and Pg R vanish in four 
dimensions, because they become full contractions of self-dual and self-antidual antisymmetric 
tensors. 

The evanescent operators listed here would look somewhat simpler if we removed from them 
all the terms proportional to e. It would be equivalent to changing one version of the MS scheme 
to another. Keeping the terms proportional to e in the above equations makes our NDR-MS 
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scheme equivalent to the one where the so-called "Greek projections" are used (see appendix 
D). 



Appendix B 

Here, we specify the evanescent operators that are necessary as counterterms for one-loop 
diagrams with insertions of the AF = 1 operators ( |3.1| ). 



£.VLL _ 




-16 + 4e)Q7 LL , 






^.VLL _ 


{r llxlvlp P L d a ){u P ^Yl P PL^) + (' 


-16 + 4e)Qy LL , 






£>VLR _ 


{s a llxlvlp p L dP){u^ri v i p p R c a ) + ( 


-4 - 4e)Q7 LR , 






^VLR _ 


(s a ^ lulp P L d a )(u^^YPR^) + ( 


-4 - 4e)Qy LR , 






^SLR _ 


(s a a^P L dP)(uPa^P R c a ) - 6eQl LR , 








£.SLR _ 


{s a a pp P L d a ){u p a^P R c p ) - QeQl LR , 








£>SLL _ 


(s^i^^PLd^YYYrPLC : 


) + (-64 + 96e)Ql LL + I 


;-16 + 1 


3e)Q 3 SLL 


^SLL _ 


(s a lplulpla P L d a ) (u l3 YYl P rP L c f3 \ 


) + (-64 + 96e)Qf L + ( 


C— 16 + i 


3e)Qf L 



The remaining evanescent operators (for the VRR, VRL, SRL and SRR sectors) are obtained 
by interchanging L and R above. 

Appendix C 

This appendix contains definitions of the "additional" evanescent operators that are not 
necessary as one-loop counterterms in the AF = 1 effective Lagrangian in section [| However, 
they have to be included before performing transformation to the "plus-minus" basis in section 

£ 3 VLL = (s Q 7 ^LC Q )(MV^L^)-Qy LL , 

£ 4 VLL = (sVLC^YP L d a )-Q2 LL , 

£ 5 VLL = {s a lplvlp P L c a ){u^^YP L d p ) + (-16 + 4e)Q^ LL , 
£ 6 VLL = (s a l^l P PLC P )(u^YYP L d a ) + (-16 + 4e)gy LL , 
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E. 
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E 



VLR 



VLR 
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3 

SLL 
4 

SLL 
5 

77.SLL 
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^SLL 



s a P R c a ){u?P L £) + )^QX LK , 
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s a a^P R c a ){u p cr^P L d p ) + 3eQ7 LR , 
S V F P/)(uV v P i( f a ) + 3eQy LR , 

^Pr^u^Pl^) + 2Qf L R , 
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^7*7^0 (^7^^) + (8 + 8e)Qp R 
^lu-ypPncP) (u^YrYPLd ) + (8 + 8e)Q 



SLR 
2 ' 



SLL 

3 > 



S a P iC «)(^P L ^) + ~Qf L - Ig 
z o 

rp^)(^p^) + iQ 2 SLL -igf L , 

ra^P L <?)(vP<T*'P L dP) - QQf LL - \qf L , 
s a a^P L c^)(u^P L d a ) - 6Ql LL - ^Qf L , 



SLL 

3 > 



^7,7p7^c*)(# 7 V7V^) - 64Qf L + (-16 + 16e)Q 
s Q 7^7p7^LC /3 )(n /3 7 ' 1 7"7 P 7 ,T J PLrf a ) - 64Q| L l + (-16 + 16e)Qf L 



The remaining evanescent operators (for the VRR, VRL, SRL and SRR sectors) are obtained 
by interchanging L and R above. 

Appendix D 

In the present appendix, the notion of "Greek projections" 0, ^, |22[ is recalled and general- 
ized to the case of SLL-sector operators. Let us denote the Dirac structure of the operator in 
eq. ( |1 . 1|) by ® Tb- The insertion of this operator in one- and two- loop diagrams results in 
new Dirac structures like 



T n T A ® T n T B , 



(D.l) 



where T n = 7 Atl 7 At2 ...7 Mn . Several examples of such structures occur in appendices A-C. It has 
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been suggested in ref. [^2j to project them onto physical operators as follows. One defines the 
projection G so that the following equality is satisfied 



G [T n T A ® T n T B ] = £G 



r 



B 



(D.2) 



In the case of Y a = Y b = 'jaPi, performing the projection G amounts to replacing £g> by 
7 T on both sides of the above equation and contracting the indices using D-dimensional Dirac 
algebra. In this manner, the coefficient £ is determined. One finds for instance: 



G 



(16 -4e) G QVll +0 ^ £ 



(D.3) 



with Q7 LL as defined in eq. (|3~l| ). 

It has been pointed out in ref. || that for a proper treatment of counterterms in two-loop 
calculations, one has to use eq. ( p.3| ) only as a prescription for defining an evanescent operator. 
In the case at hand, this is the operator E^ LL of appendix B. As discussed in ref. 0, in the 
case of VLR and SLR operators, the analogous projections are performed by replacing ® by 
1 and 7 r , respectively. Examples of the corresponding evanescent operators can be found in 
appendices A-C. 

The projections in the SLL sector are slightly more involved. In the case of the insertion of 
qSll Qr qsll^ r.h.s. of eq. ( p.2| ) has to be generalized to a linear combination of these two 
operators. The same applies to the pair (Q| LL , Qt LL )- The projection G is now performed by 
replacing (g) by 7 Q 7/3- After the projection, one finds linear combinations of g a p and "^alp on 
both sides of the equation. This allows extracting the coefficients in question. One finds for 
instance 



G 



{s a l^l P laP L d p ){u P ri v l P l a PLC a )\ = (64-96e) G fof l 1 +(16-8e) G \Qf L ]+0(e 2 ). (D.4) 



The corresponding evanescent operator is Ef hh in appendix B. An alternative approach to 



projections can be found in ref. [11 
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Appendix E 

In this appendix, the 1/e and 1/e 2 poles in the one- and two-loop diagrams are given for the 
AF = 1 calculation in the SLL sector. Analogous results for the remaining sectors can be found 
in refs. and 0. The gauge-fixing parameter A is set to unity here, i.e. the Feynman-'t Hooft 
gauge is used. 

Each insertion results in a linear combination of Qf LL , Qi LL , after subtracting the evanes- 
cent counterterms (see appendix B) or, alternatively, after performing the "Greek projections" 
(see appendix D). Table [I] gives the singularities (without colour factors) in the coefficients 
of the resulting operators, for each diagram separately. The numbering of the diagrams and 
values of the colour factors are exactly as in figs. 1, 2 and tables 1, 2 of ref. ||. The multiplicity 
factors of the diagrams are included. 

In the two-loop case, the singularities include one-loop diagrams with counterterm insertions. 
The counterterms proportional to evanescent operators are multiplied by an additional factor 
1/2, and, at the same time, the term —26c in eq. ( |2.14| ) is ignored. Correctness of such a trick 
has been justified in refs. || [10| . 

The singularities from table [I] apply for the pair (Q| LL , Q4 LL ), too. After including colour 
factors and summing the diagrams, the 1/e singularities build a 4 x 4 matrix in the basis 

{Qf LL , Ql LL , Qf LL , Q 4 SLL } 

B = ^B 1+ (^-) 2 B 2 + 0{a% (E.5) 

47T V47T/ 

from which the anomalous- dimension matrix can be obtained by means of 

7 (o)sll = _ 2 [ 2 a 1 6 ij + {B 1 ) ij \, (E.6) 
7 (i)sll = _4[2a 2< y ij - + (4)«]- (E.7) 

Here, a\ and a 2 originate from 1/e singularities in the quark field renormalization constants. 
They read 



d\ — — Cp, 0,2 — Cf 



l c ?-T N+ 2 f . 



(E.8) 



Remembering the trick applied to evanescent operators here, it is easy to verify that eqs. ( |E.6| ) 
and (|E.7| ) are equivalent to eqs. (|2.13|) and Q2.14j ) from section |. 
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Table 1: Pole parts of the one- and two-loop diagrams with insertions of Qf LL and Qf LL . The colour factors 
are omitted, whereas the multiplicity (M) is taken into account. The numbering is according to fig. 2 of ref. 
0]. While the singularities in front of the resulting Qf LL and Q2 LL are t ne same in this table, they become 
different after the inclusion of colour factors. The same comment applies to Qf LL and Q4 LL . When the colour 
factors are omitted, the results for Qf LL and Qf LL insertions are equal to those for Qf Lh and Qf LL insertions, 
respectively. 
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